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Abstract 

We study the path-integral formalism in the imaginary-time to show its validity in a case with 
a metastable ground state. The well-known method based on the bounce solution leads to the 
imaginary part of the energy even for a state that is only metastable and has a simple oscillating 
behavior instead of decaying. Although this has been argued to be the failure of the Euclidean 
formalism, we show that proper account of the global structure of the path-space leads to a valid 
expression for the energy spectrum, without the imaginary part. For this purpose we use the proper 
valley method to find a new type of instanton-like configuration, the "valley instantons" . Although 
valley instantons are not the solutions of equation of motion, they have dominant contribution to 
the functional integration. A dilute-gas approximation for the valley instantons is shown to lead to 
the energy formula. This method extends the well-known imaginary-time formalism so that it can 
take into account the global behavior of the theory. 

PACS numbers: 03.65.Db, 03.65.Sq, 11.90+t 



I. INTRODUCTION 



The imaginary-time formalism has been successful for studies of the various quantum tunneling phenomena in the 
semi-classical regime. This is because of the existence of the solution of the Euclidean equation of motion, around 
which we could evaluate the relevant functional integration in Gaussian and higher order approximations. In the 
systems with finite degrees of freedom, this is known to lead to WKB results. 

More specifically, in case when there are perturbatively degenerate vacua (or ground states) , the instanton calcula- 
tion takes into account the tunneling between them and gives the nonperturbative contribution to the energy splitting 
jl],^| . On the other hand, when a perturbative vacuum is unstable due to the tunneling to lower energy states, the 
bounce solution |^-^| leads to the imaginary part of the energy, thus the decay rate. This formalism provides us with 
a good calculational tool, valid in a wide range of physical systems. 

One subtle feature of these calculations can be seen for a one-dimensional quantum mechanical model with a 
potential, V (</>), illustrated in Figj^. (We denote the coordinate by <f> in this letter.) If one restricts the wave functions 
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to have only the outgoing component at <j> ^> 0esc> the hermiticity of the Hamiltonian is violated and the energy 
eigenvalues become complex; this imaginary part is a direct consequence of the instability of the localized wave packet 
at (j> ~ 0. In the imaginary-time formalism, this complex energy is thought to be evaluated by an analytic continuation 
of the ill-defined divergent Gaussian integral over the negative mode direction at the bounce solution |^-^| . 

We may then ask ourselves what are the eigenvalues for an unbalanced double- well potential illustrated in Fig.^J [|| . 
(In order to distinguish this case from the previous one, we call this 'metastable' case and the previous one 'unstable'.) 
As far as we stick to the classical solution and take into account only the contribution of its infinitesimal neighborhood 
for the path integral, the result is the same as the unstable case; we have a similar bounce solution with a negative 
mode and obtain complex eigenvalues. 

This imaginary part is obviously a fake. The perturbative ground state in the left well is only metastable in the 
sense that any wave packet that tunnels to the right well oscillates between the two wells. This oscillation generates 
the splitting behavior similar to the degenerate case, e = 0, but not to the decay rate. In other words, we can only take 
wave functions with decaying exponential at cfr 3> (fro for the potential in Fig.|] and the hermiticity of the Hamiltonian 
cannot be violated. 

It was claimed that Euclidean formalism is doomed due to the fact that it only uses the information around the 
infinitesimal neighborhood of the solutions As an alternative, the complex-time method was proposed, in which 
the real-time part of the classical trajectory takes into account the correct boundary condition for <fr > <fr = 4>esc- 
Although this is a very exciting development by itself ^,0], it is not clear whether it should replace the Euclidean 
formalism. 

In this letter we show that the proper treatment of the imaginary-time path integral leads us to the correct behavior 
of the energy eigenvalues for the metastable potential. The essence of the improvement is the use of the proper valley 
method, which was developed independently by Silvestrov |)| and two of the present authors [Tc| ] . It teaches us how 
to enlarge the set of background configurations besides the classical solutions in order to take into account the global 
structure of the functional space. Using this method, we construct "valley instanton" , which should replace the bounce 
solution. Interestingly, it has a zero mode and this expedites the calculation of its determinant and Jacobean as was 
the case of the instanton. It is a well-localized configuration with respect to the imaginary time and their dilute gas 
sum generates the reasonable energy corrections of the lowest states instead of the imaginary part. We will also show 
that it converges analytically to the instanton in the limit of e — > and all the results reproduce the well-known 
instanton results. 



II. THE VALLEY BOUNCE AND THE VALLEY INSTANTON 



Consider the quantum mechanical system with the following Euclidean action; 
with the potential, 

V{4>) = l -4> 2 (1 - gcfrf - e(W - 3 5 4 ^ 4 ), (1) 

where the coupling constants g and e are positive. The potential (Q) has a local minimum, V(0) = and a global 
minimum at (f> = 1/g, where V(l/g) = -e. (Fig.| is plotted for g = 0.3 and e = 0.25.) The potential V{<fr) in Eq.(@) 
is a canonical form of quartic potentials with two minima, since any such potential can be cast into this form by 
suitable linear transformations on <fr and r. In this sense the following analysis is quite a general one. In the following, 
we consider the cases with small coupling j< 1, but not necessarily small e. 

As mentioned in the introduction, the naive application of the semi-classical approximation around the bounce 
solution of this potential leads us to the fake imaginary part of the energy eigenvalues. In order to circumvent this, 
we construct the different type of configurations around which we expand the action. Such configurations are most 
straightforwardly identified by the valley methods [|| , or more specifically, by the proper valley method . The 
latter has many advantages over the previous one. It has been applied to the Borel summability problem pj, induced 
bubble nucleation problem |T^ ] , the instanton in gauge-Higgs system p3| , and the baryon number violation problem 
Jlij successfully. 

For the current model, the advantage of using the proper valley method is that we can see how the action behaves 
if we go along the negative mode direction at the bounce, the direction most important to evaluate the path integral 
correctly. As we will see below, it reaches to <fr ~ (fro, and see the behavior of the potential around there. 
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The proper valley configurations are given by the following new valley equations; 

-dfy + V'M) =F, (2) 
(-a? + V" (</>)) f = XF. (3) 

For F = 0, this set of equations reduces to the ordinary equation of motion. Thus, any solution of the equation 
of motion is the solution of the new valley equations. Otherwise, eliminating F one finds that the left-hand side 
of Eq.(^), which is 8S/8(f>, is the eigenvector of the 5 2 S/54>84>(= D) with the eigenvalue A. The general solutions 
can be parametrized by the eigenvalue A, or any arbitrary function a of A. We denote this one-parameter family 
of the solutions, the valley trajectory, by Q (r). The contribution of such a configuration to the vacuum transition 
amplitude, 



(</> = O,t = +oo|0 = O,t = -oo) =7V / X>0e _SW>] , (4) 

J0(±oo)=O 

is evaluated around this valley by inserting the following triviality; 



1 = &J daS^J dr(0 - ■ (5) 

This forces the expansion of <f> — <j> a in the subspace orthogonal to 6S/5(f>, enabling the Gaussian integration without 
the linear terms. At the leading order of h, we obtain; 

(^ = O,t = +oo|0 = O,t = -oo) = N { da - 1 J a e- S[M , (6) 

J v det D' 

where det D' is the usual determinant less the eigenvalue A in the new valley equation (||) and 



dr^F 



dS[4> 



•/, ' , : ^ = da . (7) 

drF 2 J J drF 2 

This way we are able to perform the integral for the negative mode direction on a more rigid basis than the subtle 
Gaussian integral. The factor J a is the Jacobian for this change of the integration variable. We also note that Eq. (|^) 
is apparently invariant under any local reparametrization in a. 

We have carried out the numerical analysis and obtained the solutions plotted in Fig]|. This is a simple extension of 
the previous analysis for the meta-stability problem in a quantum field theory by two of the current authors, H. A. and 
S. W. j^J. The solid line, a, is the bounce solution of the equation of motion. The rest have F ^ do not satisfy 
the equation of motion. We call these solutions (including the bounce solution) "valley bounce" . The values of the 
action, S, and the eigenvalue, A, for the valley bounces are plotted in Fig.[|. The bounce solution lies at the top of 
the line of the action in Fig.[|, corresponding to the fact that it has a negative eigenvalue. This negative eigenvalue 
can be read from the corresponding point of the plot of A in the lower half of Fig.|]. The most notable feature is that 
the large size valley bounces have a clean interior, where 4> = <^o & n d F = 0. (This property is shared by the higher 
dimensional configurations, i.e., the valley bubbles.) The effect of this is apparent in the behavior of the action in 
Fig|| the action decreases linearly with large \<f)\, which is proportional to the size of the valley bounce. Thus, the 
valley bounces provide a natural way for the evaluation of the contribution of large regions of the true minimum <j>Q . 

In the large valley bounces, c, d, and e in Fig.[| we notice that the shape of the wall, i.e., the transition region from 
4> ~ 0o to ~ 0, is almost identical to each other. That is, they overlap with each other very well when translated 
in r. Thus these large valley bounces can be approximated by simply connecting the walls by a flat region, <j) = <fio, 
at various separation. The shape of the wall can be most readily identified when the size (\4>\ = J drip) of the valley 
bounce becomes oo. In this limit, the wall is simply a localized transition from <j) = to <p = <po (or vise- versa). Such a 
solution is an analogue of the instanton (or anti-instanton) . The difference is that now it is not a solution of equation 
of motion, but is a solution of the new valley equation. This kind of configuration is called the "valley instanton" [jl3| . 
In the following, we evaluate the properties of the valley instantons, in order to approximate the large valley bounces 
by a pair of a valley instanton (cf> = — > </>o) an d a valley anti-instanton (cf> = cf>Q — > 0). 

Since we define the valley instantons at the large-size limit, \<f>\ — > oo, of the valley bounce, the plot of the eigenvalue 
in FigJ^ implies that the eigenvalue A of the valley instanton is exactly zero. This is not a trivial property. A solution 
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of equation of motion is guaranteed to have zero modes corresponding to its symmetry transformation, such as a time 
translation. Arbitrary background configurations do not have this property in general. However, we can prove the 
existence of the zero mode as follows: Take a derivative of Eq.(||) with respect to r, multiply F, and integrate over r. 
After partial integrations (which surface terms vanish), we then find that 

/OO POO 
F4>dT = I FFdT . (8) 
-oo J — oo 

The integral in the left-hand side is generally non-zero due to the boundary conditions of the valley instantons. This 
can be seen from the behavior of </> and F in the walls in FigJ|. Since the right-hand side is zero, we find that A = 0. 
(For the valley bounces, the integral in the left-hand side is zero since 4>(t) and F(t) are even functions. Therefore, 
A 7^ is allowed.) 

We have carried out numerical investigation of the valley instanton with A = and have successfully obtained the 
solutions, which have turned out to be almost identical to the wall regions in Fig.|. Although we know of no exact 
analytical expression of the valley instanton, it can be constructed analytically for small g 2 in a manner used in the 
construction of the constrained instanton |1| as well as other types of the valley instantons jl3| . Consider the valley 
instanton in t e [—T/2,T/2] (T ^> 1). We define its central coordinate to be at the origin, r = 0, by 4>(0) = l/2g. 
The naive perturbation in eg 2 yields the following perturbative valley instanton solution; 

= ^ + 3e 3 2 T^+O((e 5 2 ) 2 ), 

(9) 

F = -6eg 2 <j>l + 36eVr<^o + 0((e 5 2 ) 3 ), 
where </>q is the instanton solution for e = 0; 

4 = -zr-^ ■ (10) 

" g 1 + e- T y ' 

From Eq.(^), it is apparent that this naive perturbation is valid only in the region close to the instanton center, 
\t\ <C l/(e<? 2 ). On the other hand, in the asymptotic regions, r — > ±oo, we linearize the new valley equation and 
find the general solutions valid for |r| 3> 1. Coefficients of the general solutions are fixed by matching it with the 
inner solution Eq.(||) in the intermediate region, 1 <C |t| <C l/(eg 2 ). We have found that this procedure can be done 
consistently. The resulting asymptotic behaviors are for r — > +oo; 



1 Wl + ^rle- 



(11) 



F ~ —6ege 

V"((j)o) = 1 + 12eg 2 , and for r -► -oo; 



-(l + 3eg 2 r) e\ 
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(12) 



F ~ -6ege T . 



This way, the valley instanton is constructed in all regions of r for small e. 

The action of the valley instanton is given by S 1 = l/6g 2 + e(— T/2 + 1/2) + 0(e 2 ). This action is divided to 
the volume part and the remaining (proper) part as S 1 = —eT/2 + S 1 . From the construction above, we find that 
S 1 = l/6g 2 + e/2 + . . .. However, there is a subtlety on this point: In the following we integrate over the position 
coordinate of the instantons and anti-instantons in the dilute gas approximation. These coordinates are originally the 
valley parameters (as) of the valley bounces (and their central coordinates). The 0(e) term in S 1 depends on the 
definition of these valley parameters, since the definition of the volume (T) is affected by it. Therefore, careful study 
of the small valley bounces is needed to fix this term. This term, however, has only the nonleading contribution. 
Therefore, we will not pursue this problem any further here. 

The Jacobian for the instanton position is given by, 
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J 1 = - F J= . (13) 
drF 2 

Since the contribution to the integration is dominated by the central region, the leading term of Eq. (|l3|) for small e is 
evaluated by the use of Eq. (^|) . The result is that J 1 = l/-\/6g 2 (l + 0(eg 2 )). The first term is the instanton action for 
e = 0. Therefore, in the limit e — > 0, the Jacobian of the valley instanton reduces to that of the ordinary instanton. 

The determinant, det£>', can be calculated by extending the Coleman's method ||, in spite of the fact that the 
valley instanton is not the solution of the equation of motion. This is due to the fact that the valley instanton 
possesses the exact zero mode F(t). We define the asymptotic coefficients F± by F(t) ~ F±e TuJ±T , where for the 
sake of notation we introduced ui+ — V"(0) — 1. After some calculation, we find that the ratio of the determinants 
is given by the following; 

det'(-ff + y'V)) _ (w+ ^_ )(r/a _^ = i r F , (u) 

det{-d*+w%) Ke > 2UJ+LO-F+F- J_ 00 1 ' 

The exponential factor is the perturbative contribution to the zero energy at the true minima, </>q. The factor k is 
the 'proper' instanton contribution. From Eq.(pd|) and Eq.(|l2]), we find that k reduces to the ordinary instanton 
determinant for e — > 0. 



III. PATH INTEGRAL IN THE "DILUTE-GAS" VALLEY INSTANTON APPROXIMATION 

Combining all factors evaluated in the previous section, we find the expression of the finite time (T) vacuum 
transition amplitude to be the following; 

(0 = O,T = ?> = O,T = O) ~ , 

((f> = 0,t = T\(p = 0,t = 0) ^ 

where the amplitude (. . .)o is for the harmonic oscillator with uj + , n is the number of the valley instanton pairs, 
and the factor a is the product of the proper contributions of the action, the determinant ratio and the Jacobian; 
a = (J 1 1 \fn)e~ s . The actual integrations over the positions of the valley instantons arc in the factors /„; 

for n = 0, 

UT)={ r r fT2n f r 2 (16) 
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where zero-energy contributions of the determinants are absorbed in e by e = e — (cj_ — lj + )/2. 

The infinite series can be summed by the use of the generating function method From Eq.([l6|), we find that 
the following differential equation is satisfied by Z(T); 

Z{T)" ~ eZ(T)' - a 2 Z(T) = 0. (17) 

Also, Z(Q) = 1, and Z'(0) = 0. Therefore, we find that 

. . k+e- k - T ~ fc_e- fc + T , x 

where, 



k± = --± ] J-+a 2 . (19) 



Thus we find that the energies of the two lowest states, E± is given by 
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Furthermore, from the coefficients of the respective exponents of Eq.(|18|), , we find that \(<p = 0|_E±)| = ±k±/(k+ — fc_) 
This is the main conclusion of this letter. There appear no fake imaginary parts in the energy spectrum. We observe 
that E± are equal to the eigenvalues of the following matrix; 

x \ 

(21) 

Furthermore, the weights of the state localized in the left well agree with Eq.(|l8|): Denoting the eigenvectors of Eq.([2l|) 
V± with eigenvalues E±, 






(22) 



This allows a simple explanation of the result. The energy spectrum we obtain is the same as the two-level system 
made of the perturbative ground state at <p — and <p = <pQ, with the tunneling matrix clement a. 

Before discussing the implication of the result Eq. (po|) , we examine the validity of the dilute gas approximation we 
used. The mean size of the bounce, which is the mean distance, R, between the instanton and the anti-instanton 
located at the right of the instanton, can be obtained as the expectation value, (t2 — t\) in the amplitude, Eq.(|l5|) 
and Eq. ([l6|) . (Any other (r2 m — T2 m -i) with integer m would yield the same result for T — > oo.) Using the same 
generating function method as above, we obtain the following; 

(23) 

Similarly, the mean distance, d, between the anti-instanton and the instanton located at the right of the anti-instanton 
((73 - r 2 )) is given by, 




1 / I I 2 

. J- n/2 



,2 



2 



(24) 



For e > a, d ~ 1/f. On the other hand, the thickness of the instanton is 0(l/\/e) for e> 1, which can be seen by a 
simple scaling argument on the new valley equations, Eq.(^) and Eq.(||). This means that the dilute gas approximation 
is valid for e < 1. 

For e -c a, the energy spectrum Eq.(po|) gives the following; 

E ± = ^±±^±a-Uo( e -). (25) 



2 V", 

In the limit e — > 0, this result reduces to the well-known instanton result for the degenerate case. In addition, we 
have the average of the perturbative zero-point energies of the left- well (uj + /2) and the right- well (u>-/2 — e). Since 
the wavefunction is distributed evenly at the zeroth order of the e expansion, this is the correct formula for the two 
lowest energy eigenvalues. Since the perturbative contribution to the energy splitting, 8E = E + — E-, is of order eg 4 , 
even a purist would retain our result, SE = 2a. For aCKl, Eq.(|2C|) leads to, 



E+= 2 + T + °\e J ) ' E -=- e+ — - — + \ — r 1 2,i 1 

The lower energy, E— corresponds to the state almost localized in the right well. 



IV. CONCLUSION AND DISCUSSION 



In this letter we have applied the proper valley method to take into account the large nonperturbative configurations 
in the path integral. We have found and constructed a new type of instanton, the valley instanton, both numerically 
and analytically. We further evaluated its action, Jacobian, and determinant and found that these have a smooth 
limit to the ordinary instanton for the degenerate case. The dilute valley-instanton approximation to the path integral 
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has lead us to the energy formula Eq.(|2(j). Thus we have successfully shown that proper treatment of the imaginary 
time formalism does not lead to any contradiction and in fact yields the valid energy formula. 

Although the dilute valley instanton gas approximation fails for e > 1, this does not limit the applicability of the 
proper valley method. For such a case, the valley bounces on the background 4>q are expected to become important. 
Such valley bounces are known to exist from the analysis of jl2| and allows us to take into account the contribution 
of the configurations with (j> < 4>o- 

In view of the current development, the ordinary calculation of the imaginary part of the energy in the unstable 
cases (Fig.pl) needs to be examined under the new light. This is under way and will be reported elsewhere. 
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FIG. 1. A potential that is flat in the asymptotic direction, V(oo) = — e. 




FIG. 2. An asymmetric double-well potential, in which cj> = is only metastable. 
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FIG. 3. Valley bounce solutions (</>(t), F(t)) of the new valley equations. Center of all of the configurations are chosen to 
be at the origin, r = 0, around which the solutions are symmetric. The solid line a in the upper figure is the usual bounce 
solution, which has F(t) = 0. The other lines, 6-e, are unique to the new valley equations. 
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FIG. 4. The values of the action, S, and the eigenvalue, A, of the valley bounces. The peak of the action is given by the 
bounce solution, the line a in Fig. 3. The points corresponding to the valley bounces in Fig. 3 are plotted with circle. The solid 
lines are drawn as the guide for eyes. The valley parameter is chosen to be \<f>\ = J drfi. 
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